Characteristics of Two-Dimensional Quantum Turbulence in a Compressible Superfluid 
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Under suitable forcing a fluid exhibits turbulence, with characteristics strongly affected by the fluid's confin- 
ing geometry Here we study two-dimensional quantum turbulence in a highly oblate Bose-Einstein condensate 
in an annular trap. As a compressible quantum fluid, this system affords a rich phenomenology, allowing cou- 
pling between vortex and acoustic energy Small-scale stirring generates an experimentally observed disordered 
vortex distribution that evolves into large-scale flow in the form of a persistent current. Numerical simula- 
tion of the experiment reveals additional characteristics of two-dimensional quantum turbulence: spontaneous 
clustering of same-circulation vortices, and an incompressible energy spectrum with k~ 5/3 dependence for low 
wavenumbers k and k~ 3 dependence for high k. 



A critical distinction between hydrodynamic turbulence in 
a bulk fluid 1 1 ] and in one whose flows are restricted to two 
dimensions is that energy dissipation at small length scales 
is generally inhibited in the latter. In two-dimensional (2D) 
flows subject to small-scale forcing, energy flux is blocked 
through the small length scales and, instead, energy is trans- 
ferred towards larger scales, comprising the inverse energy 
cascade of 2D turbulence El [3). Small-scale forcing may thus 
generate large-scale flows, as seen for instance in dispersal of 
atmospheric and oceanic particulates [4], flows of soap films 
0[6l and plasmas Q, and Jupiter's Great Red Spot MM- 
However, the nature of 2D turbulence in quantum fluids is less 
clear. Progress in 2D quantum turbulence (2DQT) may of- 
fer innovative routes to understanding quantum fluid dynam- 
ics (10| [TU and aspects of the universality of 2D turbulence. 
Here we describe an experimental and numerical study of 
forced and decaying 2DQT in a dilute-gas Bose-Einstein con- 
densate (BEC). Our primary result is the first clear evidence 
that three key characteristics of 2D turbulence may also simul- 
taneously appear in systems exhibiting 2DQT: (i) emergence 
of large-scale flow from small-scale forcing, seen experimen- 
tally and numerically, (ii) numerical observation of the forma- 
tion of coherent vortex structures accompanying approximate 
enstrophy conservation [12], and (iii) numerical observation 
of an incompressible kinetic energy spectrum with &~ 5/3 de- 
pendence for low wavenumbers k and k~ 3 dependence for high 
k. Our observations are consistent with the notion that an in- 
verse energy cascade can exist in this system. 

Concepts of significance for 2D turbulence and quantum 
fluids share a common origin. Analyzing point vortex motion 
in a bounded domain, Onsager proposed in 1949 that long- 
lived vortices may develop via mergers of smaller vortices in 
turbulent flows of a 2D fluid, enabling the remaining vortices 
to move more freely and thereby maximize system entropy 
US [14]]. He also proposed that vortices in superfluids have 
quantized circulation, and implied that turbulent 2D vortex 
dynamics might be ideally studied in superfluids. However, 
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FIG. 1 . (a) Timing sequence used to study 2DQT. (b) and (c) Exper- 
imental in situ column-density images of the BEC immediately prior 
to the stir, as viewed (b) in the plane of 2D trapping and (c) along the 
z axis. Lighter grayscale shades indicate larger column densities, as 
in subsequent experimental and numerical density data, (d) Illustra- 
tion of stirring, the black arrow shows the trajectory of the harmonic 
trap center relative to the larger fluid-free region created by the laser 
barrier (e) In situ image of the BEC 10 s after stirring; vortices are 
not observable, necessitating an expansion stage to resolve them. 



experimental demonstration of 2DQT has not been reported 
and has only recently been addressed numerically fT5lfT9ll . 

To experimentally reach the 2DQT regime, we utilize op- 
tical and magnetic confinement to create highly oblate BECs 
l20l . A harmonic potential with radial (r) and axial (z) trap- 
ping frequencies (a) r /2n, oj z /2n) = (8, 90) Hz confines BECs 
of up to ~2xl0 6 87 Rb atoms having radial and axial radii of 
52yum and 5 /mi respectively. For these conditions, vortex 
bending and tilting are suppressed |2T1L leading to 2D vor- 
tex dynamics. An annular trap is created with a 23-//m radius, 
blue-detuned Gaussian laser beam directed axially through the 
trap center; the beam creates a barrier of height Uq ~ 1.5yUo, 
where jiq ~ Moj z is the BEC chemical potential in the purely 
harmonic trap. Relative to the phase transition temperature 
T c ~ 116 nK, the initial temperature is T ~ 0.9r c . 

At time t - of each experimental run, a magnetic bias 
field moves the center of the harmonic trap, but not the central 
barrier, in one complete 5.7-yum-diameter circle over 333 ms. 
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This motion induces small-scale forcing and nucleation of nu- 
merous vortices in a highly disordered distribution, which we 
identify with 2DQT much as the notion of a 'vortex tangle' 
is identified with 3D quantum turbulence QUI El [23) . Af- 
terwards, the BEC remains in the annular trap for a variable 
hold time th up to 50 s while the 2DQT decays. At th =1.17 s, 
the system temperature is reduced to ~0.6r c in order to de- 
crease rates of thermal damping and vortex- anti vortex recom- 
bination. At the end of the hold period, the central barrier is 
ramped off over 250 ms, the BEC is released from the trap, and 
ballistic expansion of the BEC enlarges the vortex cores such 
that they are resolvable by absorption imaging. Figure [TJillus- 
trates this sequence and shows images of the trapped BEC. 

Two experimental time sequences of post-stir dynamics are 
shown in Figure [2ja) and (b), emphasizing the microscopic 
variability of vortex distributions. Just after the stir (th = ms) 
a disordered vortex distribution appears. Large-scale super- 
flow is evident after th ~ 0.33 s and with increasing th, as indi- 
cated by the large fluid- free hole in the middle of the expanded 
BEC; this flow evolves into a persistent current by th ~ 8.17 s. 
An optional 3-s hold between barrier ramp-down and BEC re- 
lease gives the vortices pinned by the central barrier time to 
separate enough to determine the circulation winding number 
about the barrier; see Supplemental Material Fig. SI (27) . 
Our experiment demonstrates that under suitable conditions 
of forcing and dissipation, a highly disordered vortex distri- 
bution can evolve into a large scale flow in an annular trap. 
However, measuring kinetic energy spectra and in situ vortex 
dynamics remain forefront experimental problems, motivat- 
ing us to utilize numerical modeling and analysis for further 
characterizing 2DQT in a stirred, trapped BEC. 

BECs admit a first-principles theoretical approach that is 
numerically tractable, enabling accurate modeling [24]. The 
physical system consists of a large non-condensate component 
close to thermal equilibrium and a BEC responding both to 
external forcing and to damping by the non-condensate com- 
ponent. Numerically, we focus on the dynamics of just the 
BEC. We simulate the experimental stirring procedure using 
damped Gross-Pitaevskii theory (25). The parameters most 
readily measured are the total atom number N and tempera- 
ture T. We have developed an efficient Hartree-Fock scheme 
for determining the chemical potential fi(N, T) and reservoir 
cutoff energy e cut (N, T) in Ref. (26), and adapt the same pro- 
cedure to the present experiment, accounting for the shift in 
the trap minimum caused by the central barrier. We thus find 
the parameters needed to model the experimental system (27) . 

Figure [2jc) shows simulations that correspond to experi- 
mental observations. Here too vortices become pinned to the 
central barrier to form a persistent current; at ^=8.17 s, three 
vortices are pinned to the barrier, as indicated by the corre- 
sponding quantum phase profile (see Movie SI (27)). Ramp- 
ing off the obstacle beam in the simulation (over 250 ms as 
in the experiment) gives the column densities shown in Fig- 
ure |2jd), with density distributions more readily compared to 
(a) and (b). Development of superflow at th = 8.17 s in (c) 
leads to a large region of low density in the trap center after 
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FIG. 2. (a) and (b) 200-//m- square experimental column-density im- 
ages acquired at the hold times t h indicated above the images. Each 
BEC undergoes ballistic expansion immediately after the central bar- 
rier ramp-down in order to resolve the vortex cores. Each image is 
acquired from a separate experimental run. (c) In situ numerical data 
(9 6-//m- square images) for the hold times indicated. See also Movie 
SI (27). For each state represented in (c), ramping off the laser bar- 
rier in 250ms gives the data shown in (d). 

barrier ramp-down, as seen in (a), (b), (d). 

Analysis of our numerical simulations further characterizes 
2DQT through two distinct dynamical features of the system 
evolution, namely the development of a logarithmically bi- 
linear kinetic energy spectrum, and the formation of tightly 
bound, long-lived clusters of vortices with the same sign of 
circulation. To examine numerically the dependence of the 
kinetic energy on the wavenumber k at any instant in time, we 
use the techniques of previous studies (T6UT91I28) for extract- 
ing E l (k), the portion of a BEC's kinetic energy spectrum that 
corresponds to an incompressible superfluid component, de- 
rived by extracting the divergence-free density-weighted ve- 
locity field that embeds vorticity; the curl-free part of this field 
corresponds to sound waves and acoustic energy. 

The spectra of Fig.[3]are obtained from various times of the 
simulation and calculated using spatial grids of 181 1 2 points 
separated by £/4 = 0.1 /mi, where the £ is the healing length. 
Each curve shows the spectrum of a 2D slice through z = 0, 
although the spectra are little changed by averaging slices 
through the BEC. The ultraviolet (large k) E l (k) oc k~ 3 re- 
gion of the spectrum is a conspicuous feature once vortices 
are present. This power law is a universal property of a quan- 
tized vortex core in a compressible 2D quantum fluid, as ana- 
lyzed in Ref. (29), occurring for k > k s = The associated 
length scale ~2;r£ thus serves to distinguish between scales 
where the system's physical characteristics are dominated by 
motion of point-like vortices (k < k s ), and those where char- 
acteristics derive from the structure of individual vortex cores 
(k > k s ). The ultraviolet power law only plays a role in the 
energy spectrum through its amplitude, which is proportional 
to the total vortex number (29) . The only mechanisms that can 
appreciably change the incompressible energy for k > k s are 
creation and loss of free vortices. 
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FIG. 3. Log-log plots of E l (k) (per atom) for times 

(128,160,181,208,331) ms over which forcing occurs, plotted 
against kg with healing length g = 0.42 iim. Vertical dashed lines 
indicate k T , k F , k s , and kg, defined in the text. Red and blue lines in- 
dicate E\k) oc k~ 5/3 and k~ 3 , respectively. A spectral peak at 181 ms 
appears at k F « 2nj{\\g). Inset: log-log plot of E l (k) vs. kg (labels 
omitted) over the same domain as the main plot. From top to bottom, 
curves show E l (k) oc /r 5/3 and k~ 3 (solid lines), and E l (k) at 331ms, 
after 14 s of free decay, and for a charge-3 persistent current. 

As stirring injects kinetic energy into the system a Kol- 
mogorov &~ 5/3 power-law spectrum develops in the k < k s re- 
gion, and is determined by the vortex configuration [29]. This 
spectrum spans a decade in &-space and is established by the 
end of the stir (~331 ms) which is also when total incompress- 
ible kinetic energy is maximal. Logarithmically bilinear spec- 
tra are also obtained after ramping off the central barrier (not 
shown). Post stir, Fig. [3] (inset) indicates a slow loss of energy 
with approximate preservation of the Kolmogorov power law. 
Eventually the system spectrally condenses via generation of 
a metastable persistent current with three units of circulation. 

Three additional wavenumbers are indicated in Fig. [3] The 
cross-sectional radial thickness of the BEC approximately 
corresponds to the length scale 25/mi = 2nlkj. At high 
wavenumbers, kg = 2n/g corresponds to the scale of the heal- 
ing length g = 0.42 jjm, the approximate size of the smallest 
features (e.g. vortex cores) supported by a BEC. Finally, a 
wavenumber k F is associated with the peak in the incompress- 
ible spectrum at ~ 181ms, as we now describe. 

In the classical theory of forced 2D turbulence [3], spec- 
trally localized forcing is related to the injection rates of en- 
strophy (77) and energy (e) density as kp = yfqfe OX In 
our numerical results, the forcing that precedes the full de- 
velopment of the logarithmically bilinear spectrum is associ- 
ated with the spectral peak at 181 ms in Fig. [3] At 208 ms 
the peak has dispersed and the spectrum is already approxi- 
mately logarithmically bilinear. We estimate the location of 
kp from the numerically computed changes in total incom- 
pressible kinetic energy, AE l « 2.9 x 10~ 3 jiq-N, and enstro- 
phy, AQ « 0.963 x 10" 3 jiQ-NIg 2 , occurring during this 27-ms 
time interval. We find kp = ^ACl/AE 1 = 0.51g- 1 * 2tt/(11£), 
shown in Fig. [3] and coinciding with the spectral peak. 

The physical mechanism for injection of energy and vor- 



ticity into our BECs involves coupling between pairs of 
opposite-circulation vortices and acoustic energy. In the stir- 
ring process, density waves are first generated in the BEC, 
most prominently where the fluid density approaches zero. 
These density waves decay to vortices (see Movie SI l27l ) 
in two characteristic ways. First, a density wave can develop 
into a localized dark soliton and then decay to a vortex dipole. 
A second injection mechanism involves the decay of a density 
wave near the central barrier into a single vortex within the 
fluid and a partner antivortex pinned by the barrier. 

Empirically we find that the length scale ~27t/kp coincides 
with the separation of phase singularities created from the de- 
cay of a localized sound pulse into a vortex dipole, visible 
in Movie SI 1271 . Examining the instances of vortex dipole 
creation from sound during the stir period, we find dipole 
lengths d in the range 6.7g to llg, suggesting an injection 
of incompressible energy near a wavenumber k ~ In Id. Dur- 
ing the stir, there is one case of a vortex dipole annihilating 
irreversibly to sound at t = 0.23 s, where the dipole length 
is d ~ 6.7 g. Two transient events during the nominal con- 
stant enstrophy period discussed below correspond to dipole 
annihilation at ~ 6.7 g and ~ 8.9^. Furthermore, the super- 
fluid density modulations preceding vortex nucleation at the 
181-ms spectral peak have a length scale of approximately 
llg - 2n/kp; see Fig. S2 l27l . Taken together, these obser- 
vations indicate that forcing involves efficient energy and en- 
strophy transfer from the compressible to the incompressible 
fluid components for wavenumbers kp < k < k s . 

The conservation of enstrophy in 2D turbulence corre- 
sponds to conservation of vortex number in 2DQT |[T7l |29l . 
Between -300 ms and -600 ms after the stir begins, the 
total vortex number is nominally constant in our simula- 
tion. When vortex-antivortex annihilation occurs, the result- 
ing sound pulses are quickly (within ~20 ms) refocused by 
the inhomogeneous density, regenerating the vortices. We 
thus identify this 300-ms period of nominally constant vortex 
number with enstrophy conservation. During this period, we 
numerically observe four instances of the formation of two- 
vortex clusters (same-sign vortices). Fig. [4] shows three two- 
vortex clusters and a vortex dipole (opposite-sign vortices) 
present 30 ms after the end of stir. This dipole exists for 20 
ms, while the longest- lived of the two- vortex clusters exists 
for 630 ms. The vortices of this cluster orbit each other 15 
times, travel together halfway around the BEC, and eventu- 
ally dissociate upon colliding with a vortex dipole; see Movie 
SI l27l and Fig. [4] We see in Fig.[2ja) and (b) that there are 
large regions free of vortices, and regions where many vortices 
are densely packed, which is indicative of clustering, but we 
are not able to measure the vortex circulation directions to di- 
rectly confirm this. Clustering was, however, experimentally 
observed in Ref. |20l in the form of dipolar clusters repro- 
ducibly generated by a moving obstacle. In the present case 
of circular forcing, numerically we observe long-lived clus- 
ters emerging intermittently within otherwise irregular flows, 
in a manner suggestive of Onsager's statistical hydrodynamics 
result (l3l . 
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FIG. 4. Numerically obtained BEC column density shown for three 
hold times. In each 9 6-fim- square image, symbols indicate clusters 
of same-sign vortices, either co-rotating (red), or counter rotating 
(blue) with the stir. Each shape represents the same cluster at dif- 
ferent times. Vortex dipoles and their propagation directions are in- 
dicated by green arrows. At t h = 648 ms, one cluster remains after 
having traveled clockwise halfway around the BEC. See Fig. S3 l27l . 

Over the tens of seconds after stirring stops, free vortices 
decay either by being pinned to the obstacle beam, damping 
at the outer BEC boundary, or annihilating with free or pinned 
vortices of opposite sign. This results in the formation of a 
persistent current, reflecting the net angular momentum im- 
parted by the stirring. The development of this superflow rep- 
resents the growth of large-scale flow originating from small- 
scale forcing in an annular geometry, and serves as a check 
of our numerical procedures and interpretations. The mean 
number of vortices (pinned and free) for th = 23 s is 3.5 in 
the experiment, and 5 in the simulation. For th = 43 s, these 
values decline to 2.5 and 3 respectively. 

Previous investigations of 2DQT have centered on numeri- 
cally obtaining kinetic energy spectra, but these have been in- 
conclusive regarding the possibility of an inverse energy cas- 
cade in 2DQT, conservation of enstrophy, and correspondence 
between spectra and vortex dynamics. In our study, simula- 
tion of the experimentally realized forcing shows the devel- 
opment of an inertial range. Additionally, small-scale forc- 
ing within the trap enables the clustering of vortices; as clus- 
tering suppresses vortex-antivortex annihilation, it provides a 
mechanism to enable enstrophy conservation in a compress- 
ible superfluid. Regarding the possibility of a compressible 
superfluid supporting an inverse energy cascade, energy flux 
calculations provide the most direct route to analyzing cas- 
cades, although such an approach for a trapped BEC with an 
inhomogeneous density distribution is an open problem (30). 
However, the following observations are consistent with the 
existence of an inverse energy cascade in our system near the 
end of the stir: (i) vortex dipole recombination is suppressed 
and thus there is little dissipation over a forcing range kp to k s \ 
(ii) E l (k) oc k~ 3 for k > k s , a range that cannot support incom- 
pressible energy flux |29l ; (iii) enstrophy is nominally con- 
served, and kinetic energy spectral developments occur pri- 
marily for k <k s \ and (iv) E l (k) oc k~ 5/3 for k < k s , a signature 
of conserved energy transfer across the associated scales. 

Our observations indicate that characteristics of forced and 
decaying 2DQT in compressible quantum fluids may be anal- 
ogous to those of 2D turbulence in incompressible classical 
fluids. In particular, growth of large-scale flow, aggregation 
of vorticity, nominal enstrophy conservation, and an energy 



spectrum with k~ 513 and k~ 3 spectral features occur with suit- 
able forcing. The vortex clusters are suggestive of Onsager's 
predictions, indicating a new link between Onsager's analysis 
of 2D point- vortex dynamics [ 13 ] and the theory of 2D turbu- 
lence initiated by Kraichnan, Leith, and Batchelor [2]. Our 
observations motivate further investigations of 2DQT, with 
future work focusing on energy fluxes, the roles of dissipa- 
tion and boundary conditions, and direct experimental obser- 
vations of turbulent vortex dynamics of a 2D quantum fluid. 
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Theoretical background. The damped GPE theory used 
in our simulations can be derived rigorously for the dilute- 
gas BEC through a detailed treatment of reservoir inter- 
actions within the Wigner phase-space representation (3T), 
by neglecting thermal noise. An approximate and practical 
stochastic Gross-Piteavskii theory can be obtained l32l by: (i) 
neglecting the particle-conserving reservoir interaction pro- 
cesses (scattering terms) that are known to be small in the 
quasi-equilibrium regime; and (ii) neglecting the weak spa- 
tial and time dependence of the damping parameter. This al- 
lows the damping parameter to be calculated a priori, once 
the reservoir parameters are known. The resulting stochas- 
tic Gross-Pitaevskii equation (SGPE) has been used to study 
spontaneous vortex formation l33l during Bose condensa- 
tion, and vortex dynamics [21, 26] at high temperature. The 
SPGPE is derived by treating all atoms above an appropri- 
ately chosen energy cutoff 6 cut as thermalized (incoherent re- 
gion) with temperature T and chemical potential yu, leading 
to a grand-canonical description of the atoms below 6 cut (co- 
herent region). A dimensionless rate y = y(T,ju,e cut ) de- 
scribes Bose-enhanced collisions between thermal reservoir 
atoms and atoms in the BEC. 

However, it is not known how to extract a well-defined 
condensate orbital from the SGPE trajectories in high- 
temperature systems containing vortices. To extract physical 
information about the condensate dynamics we nelgelct the 
thermal noise to obtain the damped Gross-Piteavskii equation 
of motion. We thus arrive at a tractable, microscopically de- 
termined equation for the condensate orbital, with a priori de- 
termined reservoir parameters. After a trivial shift of energy 
by the chemical potential ji we obtain the equation of motion 
for the wavefunction i/f(r, t) 

m <W^ = _ _ (l) 
ot 

where the operator L 

Wr> = (" ^ V 2 + V(r) + g\iKr, Ol 2 ) <A(r, 0, (2) 

is the generator of GPE evolution for atoms of mass m in 
an external potential V(r). The interaction parameter is g = 



47ifl 2 a/m, for s-wave scattering length a. In most cases the 
damping parameter is small (y <?c 1), and is typically much 
smaller than any other timescales characterizing the evolution. 

In applying our approach to modeling the experiment of the 
main text, we find the self-consistent parameters [i = 34faZ>, 
£cut = 83M>, for geometric mean co = (ojjoj z ) l/3 , to de- 
scribe a system of N = 2.6 x 10 6 atoms held at temperature 
T = 0.9T C in the combined trap, giving the damping param- 
eter y = 7.96 x 10" 4 . The potential is characterized by half- 
width (To = 16.3 yum, and is well contained within the coher- 
ent region: cr «c R cut = ^2e cut /mco 2 = 73 jam, which is the 
spatial cutoff imposed by 6 cut . It thus has no other significant 
effect on the incoherent region. These parameters give an ini- 
tial state containing ~ 6 x 10 5 coherent region atoms, with the 
remainder in the incoherent region. 

There are some technical limitations of this approach. 
First, the number of atoms in the simulations is approxi- 
mately constant due to the fixed chemical potential of the 
thermal reservoir, while the BEC number in the experiment 
first grows as T/T c is reduced to ~ 0.6, then decays with a 
l/e lifetime of 24(3) s. A related issue is the slow spatial 
drift of the barrier beam of the experiment, which can 
decrease the number of vortices that can be stably pinned 
to the barrier, providing a loss mechanism for the persistent 
current at long times; this latter aspect of our work will be 
discussed in a separate publication. However, our procedure 
is suitable for simulating the BEC conditions early in the 
stir and hold process, and we do not expect detailed mod- 
eling of the evaporative cooling stage to alter our main results. 

Movie SI: Dynamics of a forced, damped BEC. A movie 
of the damped Gross-Pitaevskii equation dynamics may be 
viewed online at [URL to be provided by Physical Review] . 
The movie shows the column density and phase profile in the 
z = plane generated for a simulation of damped Gross- 
Pitaevskii equation for the experimental parameters. A num- 
ber of events are visible, notably: 

i) t = 140 ms: The first evidence of a boundary vortex ap- 
pears at the inside boundary of the BEC. This is most 
apparent in the phase profile. 
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ii) t = 160 ms: A prominent density modulation feature de- 
velops. 

iii) t - 166 ms: This time is half-way through stir sequence. 

iv) t = 190 - 210 ms: The first vortices injected into the BEC 
via sound decay are visible. 

v) t = 220 ms: Many vortices are now found in the bulk 
superfluid. 

vi) t = 230 ms: Vortex dipole annihilation. 

vii) t = 260 ms: Two long-lived vortex aggregates have 
formed (~3 o'clock and 4 o'clock positions). 

viii) t = 390 - 410 ms: A vortex approaches the inner BEC 
boundary (~4 o'clock position) and is captured by the 
central barrier, increasing the circulation pinned to the 
barrier by one unit. In the process, a sound pulse is re- 
leased into the BEC. 

ix) t = 420-450 ms: A vortex dipole forms as two vortices of 
opposite circulation approach each other (~ 7 o'clock po- 
sition), then undergoes self-annihilation followed by re- 
vival at the outer BEC boundary. 

x) t = 480 - 520 ms: A vortex dipole forms (~ 3 o'clock po- 
sition), then undergoes self-annihilation followed by re- 
vival at the outer BEC boundary. 

xi) t = 590 - 640 ms: A vortex dipole forms (~ 8 o'clock po- 
sition), then undergoes self-annihilation followed by re- 
vival at the outer BEC boundary. 

xii) t = 600 ms: Thermalization of the sound field across the 
whole system has occurred by this time. 

xiii) t = 790 ms: Vortex dipole annihilation. 

xiv) t = 1000 ms: A two- vortex cluster, a dipole, and a free 
vortex collide. Vortex exchange occurs; the collision re- 
sults in the same vortex structures emerging. 

xv) t = 1.980 s: A vortex collision results in a vortex 
being slightly tilted with respect to the z axis (~ 4 
o'clock position), showing that this system is not strictly 
two-dimensional. Note that this vortex returns to an 
orientation along the z axis within about 50 ms. 



Experimental winding number determination. Multi- 
quantum vortices are energetically unstable, thus loose 
clustering is favorable to perfect co-location of multiple 
vortices. We make use of this energetic instability in order 
to determine the size of the persistent current formed and 
its subsequent time evolution in the experiment, as shown in 
Figure SI. By enabling extra hold time between beam ramp 
down and expansion, we can count the total number of free 
and pinned vortices for any hold time. For long-enough hold 
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th = 20 s ramp + 3s ramp + 3s 
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Figure SI. Vortices observed in experimental data. Here, BECs 
are held in the annular trap for 20 s prior to the central barrier ramp 
down, trap release, expansion, and imaging. In (a) the expansion and 
imaging procedure takes place immediately after the barrier ramp. 
The dark region in the center of the BEC corresponds to multiple 
units of circulation, and is not due to the presence of the barrier itself; 
i.e., with no circulation about the barrier, a hole would not appear in 
the BEC after ballistic expansion. In (b) and (c), BECs are held for 
an additional 3 s in the trap before the expansion imaging procedure. 
This additional hold time allows the vortices that comprise the region 
of vorticity shown in (a) to dissociate and become experimentally 
observable, (d)-(f) show similar images for t h = 30 s. 



times where the number of free vortices has dropped to much 
less than one per image (on average), the observed vortices 
can be attributed to current pinned at the barrier, particularly 
if they appear clustered about the position of the barrier as 
shown in Fig. SI. We note that observations of such regular 
structures of vortices, as shown, do not always occur after 
beam ramp down, and the vortex distribution is often more 
irregular. As stated in the main text, the mean number of 
vortices experimentally observed at th = 23 s is 3.5. If instead 
we remove the central barrier at the beginning of the hold 
period, and let the system evolve in a purely harmonic trap 
for 23 s, the mean numbers of vortices observed becomes 1.2. 
The fluid circulation is thus maintained at significantly higher 




Figure S2. Numerically obtained column density 181 ms after stir- 
ring begins. The bar indicates the scale of forcing calculated from 
the energy and enstrophy injected during the stir (see main text). 40- 
//m- square field of view. 



7 



ramp + Oms I ramp + 0.23 s I ramp + 1 .23 s 
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Figure S3. Vortex clusters seen for a main sequence hold time 
of th = ms, followed by the 250-ms barrier ramp-down, and an 
additional hold time of ms, 230 ms, and 1.23 s in the harmonic 
trap (96-//m-square images). In the left-most image, five two- vortex 
clusters are seen. Only two or three of these clusters existed at t h = 
ms prior to the beam ramp (red circle, blue triangle and square of 
Fig. 4 of the main text), while the remainder formed from vortices 
that were pinned to the beam and then released into the quantum 
fluid. By 230 ms after the laser barrier has been ramped completely 
off, one of the two-vortex clusters has acquired a third vortex (red 
circle of middle image), and these three vortices orbit each other. 
This structure eventually grows to form a loose five-vortex cluster 
shown 1.23 s after the barrier ramp. This cluster persists for about 
200 ms. The remaining clusters have dissociated by this time, but 
three new two- vortex clusters have formed (rectangles). 



levels in the annular trap, justifying a description of this state 
as a persistent current. 

Density modulations during forcing. Figure S2 shows the 
stir-induced density modulations observed in the numerical 
simulations. These density modulations decay to vortices. 
As density modulations at a wavenumber k correspond to 
compressible energy at that wavenumber, we expect an influx 
of energy into the incompressible regime at a wavenumber 
that corresponds to these modulations, as discussed in the 
main text. 

Vortex clusters and dynamics. In Fig. 3 of the main text 
we show images from our simulation that contain vortex clus- 
ters that occur after the completion of the stir in the annular 
trap. Alternatively, we may examine the distribution of vortex 
clusters in the gas after the ramp down of the stirring potential 
to determine whether such clusters might occur in the exper- 
imental system subsequent to beam ramp down. Clusters do 
indeed appear in such numerical data, as Fig. S3 shows, and 
are not limited to just two vortices per cluster. 
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